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Motivation 


There  is  significant  evidence  to  show  cold  (non-equilibrium) 
plasma  discharges  have  distinct  advantages  as  combustion 
ignition  /  stabilization  sources 

■  At  high  pressures  relevant  to  applications,  cold  plasmas 
generated  by  nanosecond  pulsing  that  result  in  streamer  like 
constricted  discharges 

■  Significant  experimental  difficulty  in  probing  the  structure  and 
properties  of  streamers  (small  length  scales,  short  time  scales) 

■  High-fidelity  computational  modeling  can  play  an  important  role 
describing  physics  and  chemistry  in  these  discharges 


Cold  (non-equilibrium)  plasma  discharge  in 
plasma  parameter  space 


Thermal  plasmas  ("Hot") 

■  Most  electrical  energy  goes  into  gas  heating 
(-10,000  K) 

■  All  species  can  be  characterized  by  the  same 
temperature  (in  thermal  equilibrium) 


Non-thermal  plasmas  ("Cold") 

■  Electrical  power  is  absorbed  by  electrons 
which  in  turn  produce  radicals  and  ions. 

■  Electrons  have  high  temperature  (-10,000  K 
and  more) 

■  Ions  and  Neutrals  remain  at  lower 
temperature  (-300-1000  K) 

■  Not  in  thermal  equilibrium 
(non-equilibrium  plasma) 

From  :  NRL  plasma  formulary 
http://wwwppd.nrl.navy.mil/nrlformulary/ 


Electron  energy  loss  fraction 


Characteristic  molecular  energies  and 
electron  energy  loss  pathways 

Characteristic  energy  (eV) 
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Reduced  electric  field  E/N  is  an  important 
parameter  for  plasmas  (1  Td  =  1017  V  cm2) 

Note:  Breakdown  threshold  ~  100  Td  (e.g. 
120  Td  for  air) 


Ref :  Starikovskiy  and  Alexsandrov,  Prog.  Energy  Comb.  Sci.,  2013 
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Approach  to  sustain  non-equilibrium  at  high 
pressures  (automotive  and  aerospace  appl.) 


In  principle  can  maintain  non-equilibrium  by 
high  discharge  voltages  (i.e.  high  E/n) 

■  (Rate  of  energy  gain  by  electrons)  >  (Rate  of 
energy  loss  to  gas  heating) 

However  at  high  pressures  non-equilibrium 
discharges  are  susceptible  to  Glow-to-Arc 
Transitions  (GAT) 


100 


1000 


E/N,  Td 


■  Discharge  instabilities  cause  gas  temperature  to 
rise  rapidly 


GAT  has  time-scale  of  ~100's  ns 

Can  sustain  non-equilibrium,  by  repeated 
pulsing  on  nanosecond  time  scales 

■  First  demonstration  in  early  2000  [Kruger  et  al. 
2002] 


Ref :  Starikovskiy  and  Alexsandrov,  Prog.  Energy  Comb.  Sci.,  2013 

Kruger,  Laux,  Yu,  Packan,  Perriot,  Pure  and  Appl.  Chem.,  74,  2002,  pp.  334 


e/n  ,  Td 


Nanosecond  pulsing  produces  enhanced  tail 
in  the  electron  Energy  Dist.  Func.  (EEDF) 


Air,  p=l  Torr,  300  K,  3  cm  gap 
(E/n  =  1140  Td) 


200  400  600  800  1000 

energy  (eV) 


•  Power  budget  for  nanosecond  pulsed  discharge  is  much  lower  than  a  DC 
discharge 


From:  Macheret  et  al. ,  IEEE  Trans.  Plasma  Sci.  30,  2002,  pp.  1301. 


Computational  challenges  for  plasma  ignition 
and  flame  spread  prediction 
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Time 

10ns  100ns  lus  lOus  lOOus  1ms 


Plasma  discharge  ^  v:  ^  v/: 

Gas  dynamic  I 

Flame  spread 


■  Multiple  physical  and  chemical  processes  with  vast  disparity  in  time  scales 

■  Complex  chemistries  with  high  degree  of  uncertainty 
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Coaxial  electrode  cold  plasma  igniter  for 
automotive  combustion  applications 


Coaxial  electrode  igniter 
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1200  rpm,  A/F=1 5.1  (0=1.0),  ADV:  20  deg.BTDC,  iso-octane 


of  spark  plug 

20deg.BTDC(ADV)  19deg.BTDC _ 18deg.BTDC _ 17deg.BTDC _ 


Relative  Density  at  discharge  onset  timing 
(compared  to  25deg.C,  1  atm) 

From:  Shiraishi  et  al.  SAE  Paper  2011-01-0660 


Single  electrode  (Corona)  excitation  for  automotive 
ignition  applications 


13 


Voltage  ~100kV 


Ref:  P.  Freen,  "Radio  Frequency  Electrostatic  Ignition  System  Feasibility  Demonstration,  EISG  Final  Report, 2005 
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Nanosecond  pulsed  ignition  of  supersonic 
combustion 
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Cathode  Anode 


•  7  kV  unipolar  pulses 

•  20  ns  pulse  width 

•  50  kHz  pulse  freq. 


Pure 
Oxygen 
Test  Gas 


Oblique  Jet  Nozzle 
(2mm  Dia.)  Cathode 


Anode 


Transverse  Jet 
Nozzle  (2mm  Dia.) 


OH  PLIF 


Ref:  H.  Do,  M.  G.  Mungal  and  M.  A.  Cappelli.,  "Jet  Flame  Ignition  in  a  Supersonic  Crossflow  using  a  Pulsed 
Nonequilibrium  Plasma  Discharge,"  IEEE  Tran.  On  Plasma  Sci.  Vol.36  ,  2008,  pp.  2918-2923 
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Approach 

| - - - 

■  High  fidelity  multi-dimensional  computational  simulations  of  the 
plasma  processes  relevant  to  plasma  assisted  combustion 

■  Self-consistent  plasma 

■  Multi-species 

■  Multi-temperature 

■  Gas-phase  kinetics 

■  Surface  kinetics 

■  Plasma  model  +  Gas  dynamic  model 

■  Two-way  gas  dynamic  /  plasma  coupling 


Plasma  model 


Species  continuity 


Ideal  Gas  Law 

Drift- Diffusion 
approximation 
with  bulk  convection 


Poisson's  equation 

Electron  Energy  Equation 

—  —  —  — 

+  V  •  ((j/j  E  +  V)e  -  k  V e  ) 
dt 
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Plasma  model 


■  Gas  Energy  Equation 

■  Ions  and  Neutrals  have  temperature  Tg 

■  Tg  assumed  constant,  or  obtained  by  solving  Gas  Energy 


8  y\l1[;hk 

^  -  -  ^  2m 

JisSr - +  V  •  ( YJkh  -  2Xvrg)  =  77Th(-e£/t  .V^)  +  -Me  — 

at  keH  keH  keH  Z  mkh 


Ve  ~Tg)VkJCb -e^AEfr, 


■  If  plasma  model  is  solved  with  flow  model,  Tg  is  obtained  from 
Navier-Stokes  solver  and  only  source  terms  are  calculated  by 
Gas  Energy  module 


Flow  model  (Compressible  Navier- Stokes) 
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Integral  Model  (Zheleznyak  et  al  1982): 
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Absorption  function: 
exp 
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3-term  expansion  approach  : 

I 7%h  -  (ijPoifs^  =  -AjP202I(t) 

(j  =  1,2,3) 

Sph(f)  =  Sph.  +  Sph  +  Sph 


Aj  (cm1  Torr  _1) 

Aj  (cm1  Torr1) 

Jph 

0.0067 

0.0447 

C2 

Jph 

0.0346 

0.1121 

c3 

Jp/i 

0.3059 

0.5994 

*  Luque  A,  Ebert  U,  Montijn  C  and  Hundsdorfer  W  2007  Appl.  Phys.  Lett.  90  08150 

+  Bourdon  A,  Pasko  NP,  Liu  NY,  Celestin  S,  Segue  P  and  Maroude  E  2007  Plasma  Sources  Sci.  Technol.  16  656 


Mathematical  approach  to  coupling  plasma 
and  flow  physics 
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Electrostatic  Force  Term  (No  Magnetic  field): 
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Numerical  approach 


■  ID,  2D,  3D 

■  Fully  unstructured,  hybrid  mesh 

■  Finite-volume  spatial  discretization,  backward  Euler  time 
discretization  (formally  1st  order  in  space  and  time) 

■  Flow  model: 

■  AUSM  family  of  spatial  discretization 

(2nd  order  accuracy  through  gradient  reconstruction) 

■  4th  order  RK  time  integration 


Domain  decomposition  parallel  enabled 


Plasma  chemistry  mechanism 

■  Methane-air  plasma  chemistry  mechanism 

■  Species  and  pathways  relevant  to  plasma  time  scale  (~10's  ns) 


■  26  Species  : 

E,  0,  l\l2 , 02  ,  H  ,  N2+  ,  02+  ,  N4+  ,  04+  ,  02+N2 , 02- ,  O’  ,  02(al)  ,  02(bl)  ,  0 
N2(A)  ,  N2(B)  ,  N2C  ,  N2(al),  CH4 ,  CH3 ,  CH2  ,  CH4+  ,  CH3+  ,  CH2‘  ,  H‘ 


■  85  Reactions  : 

1)  electron  impact,  2)  ion-ion,  3)  ion-neutral,  4)  neutral-neutral 
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Methane-air  plasma  mechanism 
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Methane-air  plasma  mechanism 


cluster  ion 
formation 

charge  exchange 

attachment 
ion-ion  recomb. 


Neutral  reactions  - 


dissociative 
charge  ex. 
charge  exchange 

dissociative  recomb. 


Rxn 

Reaction 

A 

B 

c 

Activation  energy 
(eV) 

Ref. 

G44 

N2+  +  N2  +  M  ^  N4+  +  M 

5.0e-41 

0 

0 

1.0 

(24) 

g45 

N4+  +  Oz  ^  Oz+  +  2N2 

2.5e-16 

0 

0 

-3.51 

(24) 

n2+  +  o2  ^  o2+  +  n2 

1.04e-15 

-0.5 

0 

-3.51 

(24) 

G47 

Oz+  +  2N2  02+N2  +  N2 

8.1e-38 

-2.0 

0 

- 

(24) 

G48 

Oy+Uy  +  N,  0,+  +  2N7 

14.8 

-5.3 

2357 

- 

(24) 

g49 

o2+n2  +  o2  ^  o4+  +  n2 

1.0e-15 

0 

0 

- 

(24) 

G50 

02+  +  02  +  M  04+  +  M 

2.03e-34 

-3.2 

0 

- 

(24) 

G51 

E  +  202  02'  +  Oz 

6.0e-39 

-1.0 

0 

-0.43 

(24) 

G52 

Oz  +  04+  302 

1.0e-13 

0 

0 

-11.64 

(24) 

G53 

Oz  +  04+  +  M  ^  302  +  M 

3.12e-31 

-2.5 

0 

-11.64 

(24) 

g54 

02  +  Oz+  +  M  ->  202  +  M 

3.12e-31 

-2.5 

0 

-11.64 

(24) 

G55 

0-  +  o2+  ^  0  +  o2 

3.464e-12 

-0.5 

0 

-10.61 

(24) 

G56 

N2A  +  Oz  ->  N2  +  20 

1.7e-18 

0 

0 

-1.05 

(25) 

G57 

N2A+02  N2  +  02(bl) 

7.5e-19 

0 

0 

-4.54 

(25) 

G5S 

N2A+  N2(A)  N2  +  N2(B) 

7.7e-17 

0 

0 

-4.99 

(25) 

g59 

N2A  +  N2(A)  ^  N2  +  N2(C) 

1.6e-16 

0 

0 

-1.31 

(25) 

G 60 

N2(A)  +  n2  ^  n2  +  N2(B) 

1.0e-16 

0 

1500 

-0.32 

(25) 

n  1 

N2(A)  +  0  ^  n2  +  0 

3.0e-17 

0 

0 

-6.17 

(25) 

g62 

N2(B)  +  Oz  ^  N2  +  20 

3.0e-16 

0 

0 

-2.23 

(25) 

G63 

N2(B)  +  n2  ^  N2(A)  +  n2 

1.0e-17 

0 

0 

-1.18 

(25) 

G 64 

N2(al)  +  02  ^  N2  +  20 

2.8e-17 

0 

0 

-3.28 

(25) 

Ms 

N2(al)  +  N2  ->  N2  +  N2 

2.0e-19 

0 

0 

-8.4 

(25) 

G66 

N2(C)  +  02  ^  N2  +  20 

3.0e-16 

0 

0 

-5.91 

(25) 

G67 

N2(C)  +  N2  ^  N2(al)  +  N2 

1.0e-17 

0 

0 

-2.63 

(25) 

G68 

N2(C)  ->  N2(B)  +  hv  (photon) 

3.0 

0 

0 

- 

(25) 

G69 

N2(A)  +  ch4  ^  n2  +  ch4 

3.0e-21 

0 

0 

-6.17 

(25) 

Gjo 

N2(B)  +  CH4  N2(A)  +  ch4 

2.85e-16 

0 

0 

-1.08 

(25) 

Gn 

N2(B)  +  CH4  N2  +  CH3  +  H 

1.5e-17 

0 

0 

3.15 

(25) 

G72 

N2(al)  +  CH4  N2  +  CH3  +  H 

3.0e-16 

0 

0 

2.1 

(25) 

G73 

N2(C)  +  CH4  N2  +  CH3  +  H 

3.0e-16 

0 

0 

-0.8 

(25) 

G74 

02*  +  CH4  02  +  CH3  +  H 

3.0e-21 

0 

0 

- 

(25) 

G75 

02*  +  02  ^  Oz(al)  +  02 

1.86e-19 

0 

0 

-3.52 

(25) 

G76 

02*  +  02  ^  02(bl)  +  02 

8.1e-20 

0 

0 

-2.87 

(25) 

G77 

o2*  +  o2  o2  +  o2 

2.3e-20 

0 

0 

-4.5 

(25) 

Gjs 

o2*  +  0  ^  o2  +  0 

5.0e-18 

0 

0 

-4.5 

(25) 

G79 

02*  +  0  Oz(al)  +  0 

2.7e-18 

0 

0 

-3.52 

(25) 

Gso 

02*  +  O  ^  02(bl)  +  O 

1.35e-18 

0 

0 

-2.87 

(25) 

G81 

N2+  +  CH4  ^  N2  +  CH3+  +  H 

1.3e-15 

0 

0 

- 

(25) 

G82 

ch4+  +  o2  ^  ch4  +  o2+ 

5.0e-16 

0 

0 

- 

(25) 

G83 

E  +  CH4+  CH3  +  H 

2.95e-12 

-0.5 

0 

- 

(25) 

G84 

E  +  CH4+  CH2  +  2H 

2.95e-12 

-0.5 

0 

- 

(25) 

G85 

E  +  CH3+  CH2  +  H 

6.06e-12 

-0.5 

0 

- 

(25) 

Reaction  Rate  Coeff.  [m3  s'1] 


Electron  impact  reaction  rate  coefficient 
computed  using  off-line  Boltzmann  solver 

Bolsig+  (Hagelaar  and  Pitchford,  2005) 


k  =  k(  E/N) 
Te  =  Te(E/N) 


k  =  k(Te) 


(Recovers  non-local  aspects  of 
electron  energy  transport) 


—  *02rot 

—  •  02vib 

02al 

02bl 

—  -02  dissociation 

—  -02  dissociative  excitation 

02  electronic  excitation 

—  -02  dissociative  attachment 

02  ionization 
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Ref:  Hagelaar  GJ  M,  and  LM  Pitchford,  Plasma  Sources  Sci.  Technol.,  Vol.  14,  2005,  pp.  722. 
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Coaxial  electrode  Nanosecond 
Pulsed  Plasma  (NSP) 


Reference: 

D.  Breden,  L.  L.  Raja,  C.  A.  Idicheria,  P.  M.  Najt,  and  S.  Mahadevan,  "A  numerical  study  of 
high-pressure  non-equilibrium  streamers  for  combustion  ignition  application Journal  of 
Applied  Physics,  Vol.  114,  2013,  pp.  083302-1-14. 
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H 


Coaxial  electrode  NSP  discharge 


High  Voltage 
Pulse  Generator 


Central 

Electrode 


Cylindrical 
Ground  Electrode 


5 

d£ 

0) 

U 

<v 

I 


60 

50 

40 

30 

20 

10 

0 

-10 

-20 


Volt 

397  \ 

Energy 

v* 

Afi  \ 

r  nrront 

Describe  initial  plasma  kernel  formation  stage 

■  ~  10's  ns  of  physical  time 

Experimental  observations 

Unbranched  streamers  propagate  from  inner  high-voltage 
electrode  to  outer  ground  electrode 
Streamer  dia  (sub-mm) 

Brighter  discharge  near  inner  electrode 
Flame  spreads  from  inner  electrode  to  outer  ground 
electrode 


100 


400 


200  300 

Time  [nsec] 

<Streamer  Discharge> 

Ref:  Shiraishi  et  al.  J  Phys.  D:  Appl.  Phys.,  42, 
(2009)  135208. 


Ref:  D.  Singleton,  S.J.  Pendleton  and  M.  Gundersen, 
J.  Phys.  D:  Appl.  Phys.,  2011,  Vol.  44,  022001. 


Coaxial  electrode  NSP  discharge  simulation 
conditions 
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■  Simulation  conditions: 

■  10  atmospheres 

■  700  K  fixed  gas  temperature 

■  40  kV  applied  voltage  (E/n  ~  143  Td) 

■  lean  A/F  ratio  (40:1  air/methane) 


^electrode 


40  kV 


t=0 


>  t 
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Coaxial  electrode  NSP  plasma  simulation 
domain 


SYMMETRY 

BOUNDARY 


I 

r  =  4  mm 


SYMMETRY 

BOUNDARY 


PLASMA 

SUBDOMAIN 


ROUGHNESS  ELEMENT 
(POWERED  ELEC  ERODE) 


1.2  mm 


DIELECTRIC 

SUBDOMAIN 


GROUND 


60,000  cells 


Simulation  domain  :  sector  of  circle 

■  20  deg.  sector  angle 

■  Characteristic  size  for  single  streamer 
propagation 

■  Roughness  element  on  inner  electrode  to  pin 
location  of  streamer 


24  processor  partition 
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Sensitivity  to  roughness  element 
configuration 

Conditions:  P=10  atm,  Tgas=700  K,  40  kV,  40:1  A/F  ratio  (lean) 


1  roughness 
element 


8  roughness 
elements 


Verified  insensitivity  to  roughness  element  configuration 


Verified  characteristic  sector  angle  for  single  streamer 
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Time  evolution  of  electron  density  and 
temperature  for  coaxial  electrode  NSP 


Conditions:  P=10  atm,  Tgas=700  K,  40  kV,  40:1  A/F  ratio  (lean) 


Electron  density  transient 

8.5  ns 


8  ns 


9  ns 


Electron  temperature  transient 

8.5  ns 


8  ns 


9  ns 


9.5  ns 


2  ns  induction  time  (defined:  time  to  reach  threshold  of  1019  nr3) 

■  Streamers  bridge  electrode  gap  in  about  10  ns 

■  Ne(peak)~  1021  nr3 ,  Te(head)  ~4eV,  Te(body)  ~leV 

Secondary  streamer  (electron  attachment  luminosity?  Self-sustaining?) 


7  ns 


10  ns 


7  ns 


10  ns 


6.5  ns 


10.5  ns 


6.5  ns 


10.5  ns 


3  ns 


3.5  ns 


2  ns 

Electron  #  Density  [m"3] 

e+21 

le+20 
ile+19 


4  ns 


i 

le+18 


4.s  ns 


4  ns 


3  ns 


2  ns 

Electron  Temp  [K] 

40000 

I 

30000 


3.5  ns 


i20000 

i 

10000 


6  ns 


5.5  ns 


5  ns 

4.5  ns 


6  ns 


5.5  ns 


5  ns 


Reduced  electric  field  profiles  along  axis  of 
coaxial  electrode  NSP 

Conditions:  P=10  atm,  Tgas=700  K,  40  kV,  40:1  A/F  ratio  (lean) 


■  Recall  breakdown  E/n  about  120  Td  (for  air) 

■  Head  of  streamer  has  significant  over-voltages  (~  500  Td)  ->  high  Te 

■  Body  of  streamer  has  no  sustaining  E-field  (E/n  ~  10  Td)  ->  low  Te 
Secondary  streamer  formation  at  end  of  pulse  with  E/n  ~  200  Td 
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Species  yields  for  coaxial  electrode  NSP 
(volume-averaged  at  9.5  ns) 

Conditions:  P=10  atm,  Tgas=700  K,  40  kV,  40:1  A/F  ratio  (lean) 


Positive  ions 


1.2E+20 


1.0E+20 


8.0E+19 


6.0E+19 


4.0E+19 


2.0E+19 


O.OE+OO 


Negative  species 

2.5E+20  i 

2.0E+20 

1.5E+20 

1.0E+20 

5.0E+19 

0.0E+00 

^  O' 


Radicals 


Charged  species  (~1020  nr3) 
Dominant  radical  0  (~1022  nr3) 
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Time  evolution  of  radical  densities  and  for 
coaxial  electrode  NSP 

■  Conditions:  P=10  atm,  Tgas=700  K,  40  kV,  40:1  A/F  ratio  (lean) 


0  radical  density  transient 

8.5  ns 


8  ns 


9  ns 


7.5  ns 


4.5  ns 


O  Density  [m-3] 

le+23 


le+22 
le+21 

I 

le+20 


4  ns 


9.5  ns 


10.5  ns 


2  ns 


2.5  ns 


3.5  ns 


Time  [ns] 

Secondary  streamer  has  significant  impact  on 
overall  radical  yield 


0  radical  distribution  in  coaxial  electrode 
NSP  at  end  of  transient 
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Significant  non-uniformity  in  0 
radical  distribution 

■  ~1023  rrr3  at  inner  electrode 

■  Consequence  of  secondary 
streamer 


0  radical  concentration  is 
evidence  for  experimentally 
observed  flame  spread  profile  ? 


Corona  ignition  -  point  to  plane 

at  infinity 


Reference: 

D.  Breden,  L.  L.  Raja,  C.  A.  Idicheria,  P.  M.  Najt,  and  S.  Mahadevan,  "A  numerical  study  of 
high-pressure  non-equilibrium  streamers  for  combustion  ignition  application Journal  of 
Applied  Physics,  Vol.  114,  2013,  pp.  083302-1-14. 
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Corona  igniter 


SYMMETRIC 

i 

i 


4  cm  — 


GROUND 


SYMMETRIC 


< - 

(plasma) 


0-cm — 


POWERED  ELECTRODE 
115  kV 


■  Simulation  conditions: 

■  10  atmospheres 

■  700  K  fixed  gas  temperature 

■  115  kV  applied  voltage 

■  lean  A/F  ratio  (40:1  air/methane) 


V, 


electrode 


t=0 


-  115  kV 

*  t 
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,  Transient  evolution  of  electron  density 

•  Conditions:  P=10  atm,  Tgas=700  K,  115  kV,  40:1  A/F  ratio  (lean) 

1  ns  5  ns  10  ns  15  ns  20  ns  25  ns  30  ns 


Peak  electron  densities  in  streamer  head  (~1021  nr3) 
Electron  attachment  in  body 


Reduced  electric  field  profiles  along  axis  of 
coaxial  electrode  NSP 

Conditions:  P=10  atm,  Tgas=700  K,  115  kV,  40:1  A/F  ratio  (lean) 


Distance  from  Prong  fm] 

■  Recall  breakdown  E/n  about  120  Td  (for  air) 

■  Head  of  streamer  has  significant  over-voltages  (~  500  Td)  ->  high  Te 

■  Body  of  streamer  has  no  sustaining  E-field  (E/n  ~  10  Td)  ->  low  Te 
•  No  secondary  streamer  formation 


40 


41 


Species  yields  for  single  electrode  geometry 
(volume-averaged  at  30  ns) 

Conditions:  P=10  atm,  Tgas=700  K,  115  kV,  40:1  A/F  ratio  (lean) 


Positive  ions 


Negative  species 


Radicals 


1.4E+22 

1.2E+22 

1.0E+22 

8.0E+21 

6.0E+21 

4.0E+21 

2.0E+21 

0.0E+00 


<VV  +# 


Charged  species  (~1020  nr3) 
Dominant  radical  0  (~1022  nr3) 


Comparison  of  species  yields  for  Corona  and 
Coaxial  electrode  geometries 
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ru 

c 

2 

8 


<D 

~G 


fD 

8 


Positive  ions 


6.0E+19 

5.0E+19 

4.0E+19 

3.0E+19 

2.0E+19 

1.0E+19 

O.OE+OO 

C 

1.2E+20 

1.0E+20 

8.0E+19 


± 


-X  .X  _X  .X 
,<$>  +*  & 


S 


4.0E+19 

2.0E+19 

O.OE+OO 


Li 


Negative  ions 

1.2E+20 


1.0E+20 

8.0E+19 

6.0E+19 

4.0E+19 

2.0E+19 

0.0E+00 

2.5E+20 

2.0E+20 

1.5E+20 

1.0E+20 

5.0E+19 

O.OE+OO 


L 


*  o'ov 


ji 


ov  cr 


V  O'  xv' 


Radicals 


1.4E+22 

1.2E+22 

1.0E+22 

8.0E+21 

6.0E+21 

4.0E+21 

2.0E+21 

O.OE+OO 

1.2E+22 

1.0E+22 

8.0E+21 

6.0E+21 

4.0E+21 

2.0E+21 

O.OE+OO 


Vs  'V1 


■  i 


vN  *<? 


Time  evolution  of  radical  densities  and  for 
coaxial  electrode  NSP 
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Baseline  Conditions: 
0  radical  density  at  30  ns 


O  Density  [m'3  ] 
le+23 

l le+22 
le+21 

i 

le+20 


P=10  atm,  Tgas=700  K,  40  kV,  40:1  A/F  ratio  (lean) 


Corona  RF  excitation 
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Problem  statement  for  Corona  RF  excited 


plasma  igniter 

axis 


RF  excitation  : 

Freq.  ~10  MHz 
Voltage  ~100kV 


axis 


Conditions:  P=10  atm,  Tgas=700  K,  40:1  A/F  ratio  (lean), 

+90kV  ->  -80kV  ->  +80kV  pulse  train  (10  ns  each) 


Simulation  strategy  for  multi-pulse 
excitation 


46 


for  a  10  MHz  pulse 

Actual  and  Assumed  Waveforms  for  a  10  MHz  pulse 
(check  attached  spreadsheet) 


Discharge  structure  dependence  on 
excitation  polarity 

Positive  pulse  Negative  pulse 


f  *  +20kV 

.  p: 

a 

i 

1 

k  point  *  -20kV 

40  mm 

plate  h 

r 

I, ,  '  +46kV 

1  —  -  -47kv 

«  t  V  ['"•/ 

“  ©  1 

i  t  • 

r  :  %  ’  \  •' 

1 _ ..i . i..l _ 1 . . _ ...r ................................ 

i 

\ 

40  mm 

r 

f  f  1  I  I  V  1  Tf 

c 

_ 

Air  plasma 
1  bar 


■  Thin  streamers  for  positive  excitation  with  low  over-voltages 

Voluminous  glow-like  discharge  for  negative  excitation  with  low  over-voltages 

■  Streamers  for  high  over-voltages  (positive  and  negative  excitation) 

Ref:  Briels,  Kos,  Winands,  van  Veldhuizen,  Ebert,  J.  Phys  D:  Appl.  Phys.,  Vol.  41,  2008,  pp.  234004  lip. 
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Electron  density  evolution  for  excitation 
pulse  train 

Conditions:  P=10  atm,  Tgas=700  K,  40:1  A/F  ratio  (lean), 

+90kV  ->  -80kV  ->  +80kV  pulse  train  (10  ns  each) 


Electron  density 


NumberDensity_E(#/m3) 
le+12  le+16  le+20 


FlectronTemp(K) 

10300 


* 

5000  65000 


Reduced  electric  field  (E/N) 


2.7920070087  526756.4164 


0  radical  density 


6.515966e+12 


3.099112e+25 


Radical  density  evolution  at  end  of  each 
pulse 
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0  radical  density 


Axial  Distance  (m) 


03  density 


co 

O 

m- 

o 

> 

„ 

‘to  CO 

S 

o  % 

1- 

CD 

-Q 

E 


— End  of  1st  Pulse 


— End  of  2nd  Pulse 


, — End  of  3rd  Pulse 


0.002  0.003  0.004 

Axial  Distance  (m) 


0.005 


Nanosecond  pulsed  ignition  of 
supersonic  combustion 


Reference: 

D.  Breden  and  L.  L.  Raja,  "Simulations  of  nanosecond  pulsed  plasmas  in  supersonic  flows 
combustion  applications,"  A/AA  Journal,  Vol.  50,  No.  3,  Mar.  2012,  pp.  647-658. 


Nanosecond  pulsed  ignition  of  supersonic 
combustion 

I  ” 

- ►  Cathode  Anode 

•  7  kV  unipolar  pulses 

•  20  ns  pulse  width 

•  50  kHz  pulse  freq. 

Shield 


Pure 
Oxygen 
Test  Gas 


10  mm 

1 

6  mm 

1 

10  mm 

1 

Oblique  Jet  Nozzle/ 
(2mm  Dia.)  Cathode 


Anode 


Transverse  Jet 
Nozzle  (2mm  Dia.) 


Ref:  H.  Do,  M.  G.  Mungal  and  M.  A.  Cappelli.,  "Jet  Flame  Ignition  in  a  Supersonic  Crossflow  using  a  Pulsed 

Nonequilibrium  Plasma  Discharge,"  IEEE  Tran.  On  Plasma  Sci.  Vol.36 , 2008,  pp.  2918-2923 
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Chemical  reaction  mechanism 


H2-02  sub-mechanism  : 


16  Species 


e,  0+,  02+,  04+,  0-,  02-  ,  H+,H2+,  0,  H,  OH,  02,  H2,  O^D),  02(a1\)l  02 (b1Iq+) 


Assumptions: 

■  Rotational  energy  immediately  heats  bulk  gas 

■  Vibrational  energy  convected  out  of  simulation  domain 
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Geometry,  mesh,  and  operating  conditions 


Plasma  Mesh 

•  8000  cells 


0.2  mm  electrode 


Trapezoidal  Pulse 
10  ns  pulse 
2.5  ns  rise/fall  time 
6  kV  peak 


0.25  cm 


-to 


0.0002- 


0.0001- 


Mach  3 
225  Torr 
300  K 


Dielectric 


V 

A 


Powered  Electrode 


30  kOhm 


10  ns- 


>time 


Dielectric 
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Unperturbed  steady  flow 


Laminar  boundary  layer  with  lower  background  number  density 


Flat-plate  leading  edge  shock 


Electrostatic  potential  and  electron  density 
transients 


Potential 

[V] 


-1000  400  1800  3200  4600  6000 


Electrons 

[#/m3] 


1.0E+18  2.5E+18  6-3E+18  1.6E+19  4.0E+19  1.0E+20 
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Charged  and  radical  species  yields  at  end  of 
pulse 


IONS 


5.0E+20 

1.6E+20 

5.0E+19 

1.6E+19 

5.0E+18 


02+  and  02'  dominant  positive  and  negative  charge  carriers 


RADICALS 


h  2.0E+21 

■  8.0E+20 

■  3.2E+20 
W  1.3E+20 

5.0E+19 
®  2.0F+19 

H  5.0E+20 

■  1.6E+20 
||  5.0E+19 

1.6E+19 
®  5.0E+18 

h  5.0E+20 
E  1 .6E+20 
|  5.0E+19 
1.6E+19 
®  5.0E+18 


h  5.0E+20 
B  1.6E+20 
|  5.0E+19 
1.6E+19 
®  5.0E+18 


02bl 


0  dominant  radical 


Gas  dynamic  response  to  nanosecond  pulsed 
discharge 


temperature  [K] 


300  475  650  825  1000 


pressure  [Pa] 


30000  32000  34000  36000  38000  40000 
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Effort  nf  fln\A/  fiolrl  nn  Hicrhorno  rlwnami^g 


•  Lower  background  number  density  in  boundary  layer  ->  higher  E/N 

•  Confinement  of  streamer  to  within  the  boundary  layer 


100  ns 


200  ns 


Flow  carries  radicals  downstream 
over  micro/millisecond  timescales 


A  note  on  parallel  computing  for  these  class  of 
problems 
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80,000  APPJ  mesh  for  500  iterations  on  Lonestar  machine  at  Texas  Advanced 
Computing  Center  (TACC) 


Processors 


-  -  m  -  -  TOTAL  SPEEDUP 

-  -  ▲>  -  -  GAS  CHEMISTRY  SPEEDUP 

-  -  •  -  -  ELECTRON  DENSITY/TEMPERATURE  EQN  SPEEDUP 
------  POTENTIAL  EQN  SPEEDUP 


•  Problems  with  large  two-dimensional  meshes  and  large  chemistries  scales  well  to  a 
few  100  processors,  cutting  simulation  times  from  -weeks  to  -  1  day.  However 
further  improvement  in  speed  up  improvement  is  limited  by  algorithmic  bottlenecks 
(specifically  the  Poisson's  eqn). 

•  New  “parallel  friendly"  discretization  approaches  to  the  Poisson's  eqn.  are  required 
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Summary 


High  fidelity  simulations  of  cold  plasma  (streamer)  discharges  at  high 
pressure  relevant  to  real  application  are  demonstrated 

■  Self-consistent  plasma  physics,  multi-species,  multi-temperature,  gas 
chemistry,  surface  chemistry,  gas  dynamics 

■  Computationally  expensive  and  needs  large-scale  parallel  computing  to 
make  simulations  feasible 

■  Simulations  provide  insights  into  discharge  physics  and  chemistry  and 
coupling  with  gas-dynamics 

■  Extension  to  large  scale  problems  with  high-performance  computing 
requires  a  rework  of  established  computational  plasma  modeling 
approaches 


End  of  Presentation 
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Plasma  kernel  formation  with  active  radicals 
is  not  a  sufficient  condition  for  ignition 


Cold  plasma  generated  radicals  are  accompanied  with  no  additional 
gas  heating 

Do  radicals  accelerate  combustion  (chain  initiation  and  branching) 
reactions  for  ignition 

Finally  are  conditions  suitable  for  flame  spread 


Same  as  classic  ignition  kernel  problem, 
except  here  kernel  is  a  cold  radical  region, 
rather  than  hot  gas  region 


Question  :  Does  the  cold  radical  kernel 
grow  in  time  or  quench  ? 


combustible  mixture 


Preliminary  computational  modeling  of 
combustion  initiation  and  flame  spread 
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Solve  reactive  gas  dynamics  problem  assuming  an  initial  radical  lmm 
kernel 

ID  Axisymmetric  transient  problem 
1  mm  kernel  size  (~  multiple  overlapping  streamer  widths) 

No  additional  gas  heating  from  plasma 
10  atm,  1500  K,  lean  mixture  with  EGR  (A/F  20:1  +  50  %  EGR) 

1  %  of  O  radicals  (consistent  with  yield  from  streamer) 


Chemistry  Mechanism:  DRM22  with  22  species  and  105  reactions 

.  Species:  H2,  O,  02,  OH,  H20,  H02,  H202,  CH2,  CH3,  CH4,  CO,  C02,  HCO, 
CH20,  CH30,  C2H2,  C2H3,  C2H4,  C2H5,  C2H6,  N2 

■  Reactive  Flow  model: 

■  VizGIow  (without  plasma  calculations)  coupled  to  Compressible  Navier- 
Stokes  solver  (VizFlow) 


Other  approaches  may  be  considered  for 
automotive  combustion  ignition  applications 
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■  Principle  requirements  : 

■  Extended  plasma  kernel  size 

■  High  radical  yield 

■  Low  loss  (volumetric;  far  away  from  surfaces) 


Spark  HSP  RFEISorECCOS 


Sub-critical  microwave  excitation  with 
external  plasma  initiation  is  a  possibility 
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microwave 


Coax-fed  microwave  can  provide  a 
volume  filling  excitation  field 

External  plasma  initiation  can  be 
used  to  keep  microwave  E-field 
subcritical 


Microwave  excitation  concept  is  not  new  for 
automotive  ignition  applications 
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A/F 

microwave 


Igniter  erosion  concerns  with  Ikeda 
concept  can  potentially  be 
overcome  with  coax-fed  microwave 


From  :  Ikeda,  et  al,  AIAA  Paper  2009-223. 


High-fidelity  modeling  capability  available  to 
simulate  microwave  plasmas  with  VizGIow 
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Fig.  1 .  Schematics  of  the  RLSA  for  the  microwave  plasma  system. 
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Characteristics  of  large-diameter  plasma  using  a  radial-line  slot  antenna 

C.  Han,81  T.  Nozawa,  K.  Ishibasi,  H.  Kameyama,  and  T.  Morimoto 

Tokyo  Electron  Ltd..  TBS  Broadcast  Center.  3-6  Akasaka  5.  Minato-ku.  Tokyo  107-8481.  Japan 


1421  J.  Vac.  Sci.  Technol.  A  24(4),  Jul/Aug  2006  0734-2101/2006/24(41/1 421/4/523.00  @2006  American  Vacuum  Society  1421 
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Summary 


Presented  an  overview  of  non-equilibrium  plasma  physics  relevant 
to  automotive  ignition  applications 

■  Nano-second  pulsed  plasma  are  efficient  way  to  generate  non-equilibrium 
plasmas  at  high  pressures 

■  HSP,  DBD,  RFEIS  devices  leverage  this  concept  in  different  ways 

■  High-fidelity  simulation  studies  of  HSP  presented 

■  Streamers  produce  copious  amounts  of  radicals  (particularly  0  radicals) 

■  Radicals  are  concentrated  at  inner  electrode  possibly  explaining  the  dynamics  of 
flame  spread  from  these  ignition  sources 

■  Showed  initial  studies  of  long  time  scale  processes  in  ignition 

■  Plasma  radical  kernel  ->  local  combustion  initiation  ->  gas  dynamic  relaxation  -> 
flame  spread 

■  Extended  volumetric  radical  kernel  possible  with  subcritical 
microwave  +  NSP  ignition 
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Trends  in  automotive  combustion  engines 
are  driving  need  for  new  ignition  sources 

Improved  engine  efficiencies  and  stringent  emission  norms  are  driving 
new  technologies  in  automotive  combustion  devices 

Improved  efficiencies  achievable  through  1)  increased  compression  ratios 
in  1C  engines  and  2)  lean  combustion 

Lean  combustion  -> 

■  Increase  in  efficiency  (power/fuel  rate) 

■  Decrease  in  flame  temperature 

->  low  NOx 

Enabling  technologies 

■  Direct  injection  (no  air  intake  throttling  losses)  ->  just  in  time  combustion 

■  Lean  with  Exhaust  Gas  Recirculation  (EGR)  ->  low  flame  temp  ->  lower  NOx 

Technological  challenges 

■  Lean  combustion  (with  EGR)  ->  ignitability  issue  is  key  problem 
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Starikovskiy  and  Alexsandrov,  Prog.  Energy  Comb.  Sci.,  2013 


Conventional  spark  plug  based  1C  engine 
ignition 
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Combustion  ignition  via  highly 
constricted/localized  spark 

Spark  is  a  thermal  plasma  with  very  high 
sensible  temperatures  (~  1000's  K) 

-  lifetime/reliability 

Chemical  initiators  for  combustion  not 
the  same  as  in  a  cold  plasma 

Limited  control  on  plasma  yield 


Nanosecond  pulsed  and  Dielectric  Barrier 
plasma-based  ignition 


Dielectric  barrier 


Cylindrical 
Ground  Electrode 


1200 


rpm,  A/F=1 5.1  (0=1.0),  ADV:  20  deg.BTDC,  iso-octane 
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Side  electrode 
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Shiraishi  and  Urushira, 
SAE_2011-01-0660 


Variety  of  plasma  actuator  concepts  exist  for 
volumetric  and  surface  flow  control 


Flow 


i 

Actuator  location  reference 


Exposed  electrode  edge 


Corke  et  al.,  Ann.  Rev.  FI.  Mech.  (2010) 


Magnetic  Field  B 


Kalra  et  al.,  Expt.  Fluids,  (2011) 


Shin  et  al.,  AIAAJ.  (2007) 


Computational  issues  in  the  modeling  of  air 
plasma  interactions  with  flows 


Extremely  high  degree  of  time  disparity  in  component  physics 
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Spatial  stiffness  due  to  discharge  structure 


Sheaths 


Electronegative  plasma 


density  density 


Pressure  ■  1400  kPa 
Peak  =  dGkV 


Pressure  =  1 400  k  Pa 
PfeakVdtage-4 5kV 


Presstre  =  1(H.3kPs 
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Maximum  voltage  line 
of  streamer  discharge 


From:  Shiraishi  et  al.f  J .  Phys.  D.f  42,  2009,  135208. 


Photoionization  (3-term  Helmholtz  equation 
model) 


Integral  Model  (Zheleznyak  et  al  1982): 


Sph(f)  =  Iff 


Hr')g(R) 

4nR2 


dV 


i(f)  = 


P  +  Pr 


■W(f) 


g(R) 


02 


6Xp  Xmin^O —  QXp  Xmax^O 2^ 
P()2  R  / Xmin) 


Luque  et  al*  proposed  approximating  g(R)/P02  using  two  exponentials  functions  and 
expanded  by  Bourdon  et  al+  to  three  terms 


sph(f)  ~  Svh  +  Svh  +  *£ 


ph 


’ph  '  '’ph 


4nR 


AjPo2exP  XjP°2R 


The  integrals  are  solutions  to  three 
Helmholtz  equations: 


v%h  -  (*}Po2)2Sjph  =  - AjP202l(f ) 

(j  =  1,2,3) 


■ 

Aj  (cnv1  Torr _1) 

Aj  (cm1  Torr1) 

Sph 

0.0067 

0.0447 

c2 
^ ph 

0.0346 

0.1121 

C3 

Jph 

0.3059 

0.5994 

*  Luque  A,  Ebert  U,  Montijn  C  and  Hundsdorfer  W  2007  Appl.  Phys.  Lett.  90  08150 

+  Bourdon  A,  Pasko  NP,  Liu  NY,  Celestin  S,  Segue  P  and  Maroude  E  2007  Plasma  Sources  Sci.  Technol.  16  656 
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Plasma  chemistry  mechanism  used  in  studies 

Plasma  Chemistry  mechanism  relevant  to  plasma  time  scale  (~10's  ns) 


Methane-air  mixtures 

■  26  Species  : 

E,  0,  N2  ,  O2  ,  H  ,  N2+  ,  02+  ,  N4+  ,  04+  , 
02+N2,  Of,  0- ,  02(al)  ,  02(bl)  ,  02*  ,  N2(A) 
,  N2(B)  ,  N2C  ,  N2(al),  CH4 ,  CH3  ,  CH2  ,  CH4+  , 
CH3+  ,  CH2-  ,  H- 


■  85  Reactions  : 

1)  electron  impact,  2)  ion-ion,  3)  electron 
neutral,  4)  neutral-neutral 


Methane-air  with  EGR  mixtures 

■  39  Species  : 

E,  0,  N3 , 03  ,  H  ,  N3+  ,  03+  ,  N4+  ,  04+  , 
02+N2, 02‘,  O’  ,  02(al)  ,  02(bl)  ,  02*  , 
N2(A)  ,  N2(B)  ,  N2C  ,  N2(al),  CH4  ,  CH3  ,  CH2 
,  CH4+  ,  CH3+  ,  CH2-  ,  H- , 

H20,  H20+,  H2,  H+,  Hf,  OH,  0H+,  OH',  0+, 
C02,  C02+,  CO'  ,  03 

■  110  Reactions  : 

1)  electron  impact,  2)  ion-ion,  3)  electron 

neutral,  4)  neutral-neutral 

Additional  :  C02,  H20  and  03  reactions 


Comparison  of  baseline  and  With  EGR  cases 
for  HSP  discharge  streamer 


82 


Baseline  (lean  A/F  =  40:1) 


9.7  ns 


With  EGR  (A/F  =  20:1  +  50%  exhaust) 


5.7  ns 


Propagation  speeds  higher  with  EGR 


Electron  density  slightly  higher  with  EGR 
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Radical  densities  for  baseline  and 
cases  for  HSP  discharge  streamer 


^BN  umbe  rDe  ns  Ity. 
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LB  N  umbe  rD  e  ns  Ity. 
LB  N  umbe  rDe  ns  Ity. 
I  I N  umbe  rDe  ns  Ity. 


_CH20f/m33_overage 
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_H  (#/  m3)_a  ve  rage 
_0(#/  ve  rag  e 

_Q2al  (#/'m3  ^average 
_03(#/m3)_ave  rage 
.OH  (#/'  mS)_ave  rage 


With  EGR 


Baseline  With  EGR 


No  significant  changes  in  radical  densities 
for  case  with  EGR 


Case  1:  Pulse  train  of  -90kV  -> 
+90  kV  -90  kV  (gas 
temperature  700K) 

Pulse  Durations: 

1st  pulse:  7  ns 
2nd  pulse:  7  ns 
3rd  pulse:  7  ns 
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tvoiuuon  or  NumDer  uensity  or 

Electrons 


NumberDensity_E(#/m3) 
le+12  le+16  le+20 


1000000000 


5.00479e+21 
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tvoiution  or  iMumDer  uensity  or  u  <* 

radicals 


NumberDensity_0(#/m3) 

le+16  le+20 

j  i  mill  i  


le+24 


4.67761e+12 


3.1551e+25 
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tvoiution  or  tiearon 
Temperature  (K) 


ElectronTemp(K) 

10000 


5000  69536.2 
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tvoiuuon  or  Keaucea  tiectric  * 

Field  (E/N) 


Reduced_E_Field_E/N(Tcl)  Magnitude 

10  100  1000  10000 

^^lll|  l  ^  I  Ml  i  i  ii  i  i^^^^ 

2.948873819  11278.362456 
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0  Radical  Number  Density  at  End  of  Different  Pulses 


Axial  Distance  (m) 
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03  Radical  Number  Density  at  End  of  Different  Pulses 
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Voltage  Amplitude 
Comparison 


Voltage:  Streamer  Propagation 


O*  number  density  [#/m3] 

1.0E+19  2.5E+19  6.3E+19  1.6E+20  4.0E+20  1.0E+21 
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8000  V 


Higher  voltages  result 
in  stronger  Electric 
Field 


Streamers  propagate 
further  as  voltage 
increases 


Voltage :  Thermal  Effects 
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Stronger  Electric  fields  result  in  greater  ion  Joule  heating 
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Conclusions 


•O  radicals  dominant  species  in  plasma  (~0.5%  peak  mole  fraction) 

•Ion  Joule  heating  dominates  gas  temperature  increase  and  results  in 
blast  waves 

•Increasing  Voltage  increases  peak  densities,  gas  heating  and 
volume  of  plasma  formed 

•Chemistry  (electropositive  vs  electronegative  plasma)  affects  - 

■  Streamer  propagation  distance/speed 

■  Region  of  plasma  formation  (inside/outside  boundary  layer) 

■  Intensity  of  gas  heating  for  different  polarities 

•Anodic  pulses  appear  more  efficient  for  supersonic  combustion 

■  Radicals  produced  over  greater  volume 

■  Less  power  lost  to  heat  (for  02-H2) 


